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We propose a graphene hyperlens for the terahertz (THz) range. We employ and numerically 
examine a structured graphene-dielectric multilayered stack that is an analogue of a metallic wire 
medium. As an example of the graphene hyperlens in action we demonstrate an imaging of two point 
sources separated with distance Ao/5. An advantage of such a hyperlens as compared to a metallic 
one is the tunability of its properties by changing the chemical potential of graphene. We also 
propose a method to retrieve the hyperbolic dispersion, check the effective medium approximation 
and retrieve the effective permittivity tensor. 



Rapidly developing terahertz (THz) science and tech- 
nology acquired a great deal of attention in recent years 
due to an enormous potential for spectroscopy, commu- 
nication, defense and biomedical imaging 1 3 . The nat- 
ural diffraction limit, however, restricts the resolution 
of the standard THz imaging systems to about a wave- 
length, which is relatively large (300 jam in the free space 
at 1 THz). To overcome this restriction one can use 
the scanning near-field THz microscopy that allows for 
even submicrometer resolution in the scattering (aper- 
tureless) configuration 4 , but such technique is slow. An- 
other solution is to use a metamaterial lens with arti- 
ficially engineered properties, for example, a negative 
index lens 5 or hyperbolic-dispersion lens (hyperlens) 6 . 
While the negative-index material lens is far from being 
implemented into the imaging systems due to high losses 
and narrow resonant frequency range, the hyperlens has 
been experimentally demonstrated in the microwave 7 and 
optical 8 regimes. The hyperlens is able to convert evanes- 
cent waves into propagating ones and to magnify a sub- 
wavelength image so that it can be captured by a stan- 
dard imaging system, a microscope, for example. 

A hyperlens usually consists of metal-dielectric layers 
(in ultraviolet and optical ranges) or of metallic wires (in- 
frared and microwave ranges). Due to the employment 
of metal the properties of the hyperlens cannot be tuned 
after fabrication. In contrast to metal, graphene, a two- 
dime ntional material with striking electronic, mechanical 
and optical properties- 9 , supports surface plasmon polari- 
tons in the THz range 10,11 that are widely tunable by the 
change of graphene's electrochemical potential via chemi- 
cal doping, magnetic field or electrostatic gating 12 . Many 
plasmonic effects and photonic applications of graphene 
have been proposed^"— . Nevertheless, up to our knowl- 
edge, no graphene based hyperlens for the THz range has 
been proposed so far (however, there has been recently 
reported a graphene and boron nitride hyperlens for the 
ultraviolet 19 ). 

In this letter we propose to use structured graphene 
for the creation of a hyperlens in the THz range. To sup- 
port our proposal we investigate the effective properties 



of the hyperbolic graphene wire medium and then con- 
struct a hyperlens out of it. We check numerically the 
performance of a full-size three-dimensional (3D) and its 
homogenized two-dimensional (2D) analogue and demon- 
strate that it has the desired subwavelength resolution 
and magnification. 

Several requirements have to be satisfied for construct- 
ing the hyperlen s 20 ! 21 . First of all, an indefinite ma- 
terial (the permittivity tensor components have oppo- 
site signs) with strong cylindrical anisotropy should be 
used. Namely the radial permittivity e r should be nega- 
tive (e r < 0) while azimuthal permittivity sq should be 
positive (so > 0). In this case the in-plane isofrequency 
contour is hyperbolic: 



where q = k r /ko is the normalized radial wavevector com- 
ponent, ft = ke/ko is the normalized azimuthal wavevec- 
tor and ko = 27r/Ao = uo/c is the wavenumber in vacuum. 
Then the waves with ft > 1, that are evanescent in vac- 
uum, can propagate in the hyperbolic medium. Mathe- 
matically this means that for every ft there exists a real- 
valued q. Moreover, the dependence g(ft) should be as 
fiat as possible. That ensures the same phase velocities 
for all spatial components (various ft). There are two 
possibilities for satisfying this requirement: to select a 
material with either a large negative e r or a small posi- 
tive Eq. 

For high transmission propagation losses characterized 
by (Im(q)) should be as small as possible. For the waves 
with ft <C 1 the radial wavevector reduces to q « y/Eo, so 
it is primarily sq that determines losses. The incoupling 
and outcoupling of the waves to the ambient medium 
should also be efficient. For normally incident waves from 
a dielectric with a refractive index n onto the flat inter- 
face with a hyperbolic medium, the reflection coefficient 
is R = = n . ^S , so in order to minimize reflection 

n+q n+y?^' 

one has to match the azimuthal permittivity e$ with the 
permittivity e = n 2 of the surrounding medium. This 
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requirement limits the range of eq. In addition to that 
to maximize the hyperlens transmission the Fabry-Perot 
resonance condition should be satisfied 
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where R\ and i?2 are the inner and outer hyperlens ra- 
dius, respectively, A e ff is the effective wavelength and m 
is an integer number. The ratio of the radii M = R2/R1 
determines the hyperlens magnification. 

Finally, since no natural electromagnetic materials 
with a strong cylindrical anisotropy exist, artificial effec- 
tively homogenous metamaterials have to be used. That 
means that its lateral geometrical period P should be 
much (at least 5-10 times) smaller than the period of 
the wave with the highest n = n max . So, for example, 
if we wish to construct the hyperlens for the free-space 
wavelength Ao = 50/im that supports the wave with the 
highest Kmax — 5, then the lateral period of the metama- 
terial should not be larger than P max = To^~~~ = lM m - 

First we analyzed the properties of the graphene wire 
medium itself. Its unit cell is a rectangular block of di- 
electric (sd = 2.34 corresponding to the low- loss polymer 
TOPAS) of the size a x x a y x a z = 0.2 x 0.05 x 1/xm 3 
(a Xl a y <C Pmax) with an embedded graphene stripe of 
the width w depicted in FigJIJi. We described graphene 
for the simulations in CST 22 as a layer of thickness A = 1 
nm with the permittivity sq = £d + i 7 wnere a s is 
the surface conductivity of graphene^ 
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FIG. 1. (a) The unit cell of the graphene wire medium con- 
sists of a graphene stripe of the width w embedded into a 
dielectric, (b) To characterize the hyperbolic medium we cal- 
culated complex reflection R and transmission T coefficients 
for various angles of incidence ip. A block of the hyperbolic 
medium is placed between high-index n s dielectric layers. 

In order to retrieve the dispersion relation q(n) we 
simulated complex reflection R and transmission T co- 
efficients for various angles of incidence ip on a hyper- 
bolic medium slab (see Fig. [TJd) with the periodic (unit 
cell) boundary conditions. We considered TM polarized 
waves (magnetic field along the ?/-axis). The surrounding 
medium was a high refractive index ns dielectric. Then 
for each k and frequency u we can restore 



where m is an integer number. Since we work in the long 
wavelength limit, the challenging choice of the branch m 
is not an issue, it should be simply m = 0. The choice of 
the sign should satisfy the passivity condition Im(q) > 
0. Knowing the dispersion dependence q(uj, k) we can 
restore the components of the permittivity tensor e r and 
so through the linear regression analysis of the dispersion 
equation (pQ) 
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The statistical coefficient of determination R sq con- 
firms (if R sq close to 1) the linear regression q 2 (n 2 ) and 
the homogenous approximation validity. For the inves- 
tigated graphene wire medium we observed R S(l > 0.95. 
We should also note that this retrieval method is ap- 
plicable not only to the hyperbolic medium, but to any 
metamaterial and that by selecting another polarization 
and/or wave propagation direction it is possible to restore 
the whole permittivity tensor. 

An example of the restoration for the graphene stripe 
of width w = 80 nm is shown in Fig. [2j The color 
contour graphs (Fig. [2^i,b) show that q(u),K) is flat at 
low frequencies, but exhibits a resonance around 17 THz. 
Detailed investigation of the electromagnetic field behav- 
ior revealed a surface plasmon resonance of the graphene 
stripe at this frequency. The q{n) isofrequency contours 
(Fig. [2k, d) are flatter and the losses are smaller at lower 
frequencies. Finally, the radial permittivity e r has the 
Drude-like dependence (Fig. [2^) with large negative val- 
ues at the low frequencies, while azimuthal sq is positive 
and has small Im(so) (Fig- EF). Thus it is advantageous 
to select a low operation frequency for the hyperlens. 

In order to select the optimal geometrical design, we 
investigated the dependence of the wire medium prop- 
erties on the stripe width (see Fig. [3]) starting from no 
graphene (w = 0) to a full graphene coverage (w = 200 
nm). As expected, in the absence of graphene we re- 
store a constant refractive index ud = 1-53 (Fig. [3^l) 
with no losses (Fig. (3Jd) and permittivities e r = £q = n 2 D 
(Fig. [3b, d), while for the full graphene coverage a typical 
Drude metal-like behavior is observed for permittivities 
e r = eq. Changing the width from w = 80 nm, which we 
discussed above, to w = 160 nm we observe larger val- 
ues of Re(q) for the normal propagation k = (see Fig. 
[3^i) (and consequently worse coupling efficiency), larger 
losses and red shift of the resonance to / = 13 THz (Fig. 
[3)3) and larger negative permittivity e r (Fig. [3t). After 
examining several widths we selected for the hyperlens 
demonstration the width w — 40 nm (not shown in Fig. 
[3]) and the frequency 6 THz. 

To check the suitability of the effective medium ap- 
proach we simulated in the CST (time domain) the full- 
size 3D hyperlens made of graphene stripes embedded 
into dielectric (jid = 1.53). One layer of the hyperlens is 
shown in the Fig. UK a). The input and output periods, 
widths and radii were chosen as P in = 200 nm, P out = 600 
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(a) Re(q) = Re(k r )/k (b) Im(q) = Im(k r )/k 
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FIG. 2. Contour plots of restored (a) Re[q(uj,K)] and (b) 
Im[q(uj, k)] for a graphene stripe of w = 80nm show absence 
of resonances at low frequencies, but a resonance at / = 17 
THz. Looking at the q(K) for certain frequencies (5, 10, 15 
THz) we observe the flat Re\q{n)\ dependence (c) and smaller 
losses (d) for the frequency. Effective radial permittivity e r (e) 
shows the Drude-like behavior, while azimuthal permittivity 
sq (f) is positive with small losses. 



nm, Wi n = 40 nm, W ou t = 120 nm, R{ n = 15.12/im and 
^out = 45.36/im, respectively. The radii are selected to 
satisfy the Fabry-Perot resonant condition (|2]). The lay- 
ers of structured graphene are assumed to be periodic in 
the direction perpendicular to the image plane (period 
a y — 50 nm). We should note that the specified sizes are 
realistic for fabrication. Multiple graphene layers sepa- 
rated with a dielectric can be made up to the size of 30 
inches — . Structuring of multiple graphene-dielectric lay- 
ers structure can be done with focused ion beam milling 
or electron beam lithography. 

Now we show the hyperlens in action when being ex- 
cited with two sources (line magnetic currents) in vacuum 
separated with distance S = An/5 = 10/mi (see the artis- 
tic 3D view of the hyperlens in work in Fig. Hfc). In the 
presence of the hyperlens two sources are well resolved at 
the output interface as two peaks separated with 30/im 
(Fig. Hk) delivering the magnification M = R2/R1 = 3 26 , 
while in case of the homogenous dielectric cylinder (no 
graphene wires) we observe a single spot (Fig. |4]i). 

Then we compared the CST results with an equiva- 
lent 2D hyperlens simulation in COMSOL 27 (scattering 
boundary conditions) with homogenized permittivities 
s r = -20.1 + 8.5z, e e = 2.73 + 0.0029*. The COMSOL 
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FIG. 3. Comparison of the properties of the graphene wire 
medium for various stripe widths (0, 80, 160 and 200 nm). 
The radial wavevector shows larger values of (a) Re(q) (that 
means worse coupling to the surrounding medium) and (b) 
Im(q) (larger losses) for w — 160 nm compared to the width 
w — 80 nm. Also, a "more metallic" Drude behaviour of e r 
(c) and higher azimuthal permittivity ee (d) nm is observed 
for w — 160. The absence of graphene (w = 0) and full 
coverage (w = 200 nm) show fully dielectric and Drude-like 
behaviours, respectively. 



results with (Fig. 0^) and without the hyperlens (Fig. 
[4f) are in a good agreement with the CST results. A 
comparison between them is shown in Fig. 0^ where the 
wave intensity at the output interface of the lens is pre- 
sented. The intensity of the peaks in the presence of the 
hyperlens is larger than in its absence, due to redistribu- 
tion of the power. The intensity simulated with the CST 
is smaller compared to COMSOL that is caused by the 
coarser spatial discretization of the tapered wires with a 
staircase numerical mesh in the CST. In both types of 
simulations the peaks are well resolved according to the 
Rayleigh criterion. The 2D COMSOL simulation, how- 
ever, took several minutes versus the 3-days long 3D CST 
modeling. 

By making a hyperlens with larger radius R2 one can 
achieve a larger magnification. For example, selecting 
R2 = lOi^i gives the magnification M — 10, so two point 
sources with separation S = 10/im are imaged to 100/im 
(see Fig. [5^) and then can be resolved with a conven- 
tional THz camera. 

It is important to test the device performance under 
the pulse excitation. In the conventional THz time do- 
main spectroscopy setup 2 (THz-TDS), a very short (sin- 
gle cycle or even shorter) THz pulse is generated. Ex- 
perimentally testing the hyperlens in the real THz-TDS 
would mean shining the short (and therefore broadband 
in frequency) transient pulse and then scanning with 
the THz near-field microscope and collecting the time- 
dependent signal at the output. We did a similar sim- 
ulation in CST, exciting 2 sources with the Gaussian 
pulse (central frequency f c = 6 THz, FWHM = 12 THz), 
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FIG. 4. (a) A single structured graphene layer that consti- 
tutes the hyperlens. (b) An artistic 3D view of the hyperlens 
in action: the image of two line sources is magnified with 
the hyperlens. Full size 3D CST simulation of the hyperlens 
(c) in action and no hyperlens (d), when two magnetic line 
sources separated by Ao/5 = 10/xm are emitting TM polarized 
waves (magnetic field perpendicular to the plane of image). 
The CST results are in a good agreement with equivalent 2D 
COMSOL simulations of the hyperlens (e) and no hyperlens 
(f). Comparison of the intensity distribution at the output 
interface of the hyperlens (g) confirms that the images are 
well resolved. 



recording the field with the time monitor and then imag- 
ing the maximal field in each point during the simulation 
time (Fig. Eb). Since the graphene hyperlens is not based 
on a resonant medium, it can operate in an extended 
range of frequencies and two sources are well magnified 
and resolved (Fig. [5b). 

Due to reciprocity the hyperlens can be used not only 
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for imaging, but also for THz power concentration into 
a small volume. We wish to note that the employment 
of metal for the considered hyperlens design is hardly 
possible. In order to obtain the same conductivity of the 
unit cell as of the regarded graphene stripes, the cross- 
section of the metallic wire (for example, silver- 8 ) has to 
be of 2/im 2 . Fabricating and arranging so thin and long 
metallic wires into the required pattern is beyond the 
possibilities of the current nanofabrication technologies. 




FIG. 5. (a) A thicker hyperlens with R2 = 10Ri magnifies 
two subwavelength sources until the images can be captured 
with a conventional THz imaging system, (b) Tracing the 
maximum of the broadband THz transient pulse shows that 
the hyperlens works in the range of frequencies around 6 THz. 

We should emphasize that the scaling up the metallic 
wires together with the unit cell is not possible, since the 
period of the hyperlens should be subwavelength even for 
the higher-order spatial harmonics. Another important 
advantage of the graphene hyperlens compared to the 
metal based one is its tunability by the graphene chemical 
potential change. Thus it is possible to make the device 
reconfigurable and to resolve subwavelength features or 
concentrate THz pulses on demand. 

In conclusion, we have shown that structured graphene 
layers embedded into dielectric (graphene wire medium) 
can be used to create a hyperlens. We have proposed 
the realistic geometrical design for the hyperlens for the 
THz radiation and proved that it can resolve two line 
sources separated by a distance Ao/5. We also showed 
that time-consuming 3D simulations are in a good agree- 
ment with the quick homogenized 2D hyperlens model- 
ing, which simplifies the hyperlens engineering. 
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